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for MHz Band Wireless Power Transfer System 





IAI ee 


Keita Asano, Toru Uno, Takuji Arima 


Tokyo University of Agriculture and Technology 


Outline 


|. Introduction 
2. Proposed method 

2.1 ARMA/FDTD Algorithm 

2.2 Estimation method of Electromagnetic field 

using a transfer function 

3. Simulation result 

3.1 Calculation of the steady state E field 

3.2 Calculation of the Reflection coefficient 


4. Conclusion 





1. Introduction 
FDTD method 


Electromagnetic field analysis method in the time domain 


Advantage 
Easy modeling , Good PC memory efficiency, 
suitable for vectorization and parallelization 


Jb FDTD Application 


EM field analysis of the Products of complex structure and Human body 
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Products of complex structure 
(Smartphone, TabletPC, etc) Human body and the antenna 
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"FDTD Disadvantage 
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Human body — WPT system(kHz or MHz band) 


Need a small cell for modeling Low-frequency EM field Application 


3 By the constraint that Cournt's stability condition 


__Ttrequires Long Calculation time | 
es | A 


For example, Calculation of low-frequency EM field problem 
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Analysis condition i6 
Cell size : 5mm cubic cell 
Ti : dt-9.5 
ime step : dt-9.5ps Calculated 
en By FDTD time[us] 
Low frequency wave , 
(MHz or kHz band sin-wave) Time until steady state : T=18ps 

Time step : dt=9.5ps 


The number of time steps 


= T/dt 
1,900,000 steps 





Voltage| V ] 
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tme[us] 
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Purpose of This study 


Electromagnetic field analysis in low frequency 
Analysis of the antenna characteristics 
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Two proposed method 
1. Estimation of Rational functions with the Auto-Regressive Moving Average(ARMA) 


fi 





Value of the time domain L— —» Š 3 
° = =] 
x(#) : Input signal (OUEST à | or &| 
y(£) : Output signal by ARMA ` | N | 
Frequency|Hz| Frequency|Hz] 
Input Impedance Transfer function 


2. Estimation method of Electromagnetic(EM ) field using a transfer function 
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pulse Transfer function 


Calculation contents in this study 





Feed Point Resistor =50 Q 


WPT system near layered dielectrics 


Tested for two types of problem 


1. Steady state E field in time domain 


2. Reflection coefficient of loop antenna in frequency domain 


iL We will evaluate 


Calculation accuracy and The number of time steps 1 
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2. Proposed method 
2.1 ARMA/FDTD Algorithm 


To approximate the rational function in the frequency domain 


Y(Z) ao-taZ '-taj*t---ajZ^? 


H(Z) 2 —— = —— a t 
2 X(Z) 1-c*hZ!-ctb*--:-tbjzZ? 


p,q : The order of the rational function model 
T : Sampling time 


$ Solve the (2) in order to determine 


Equation of ARMA model 
D q 
yG) 2 — » biyn - D » ajix(n-j) 2 
i=1 j=0 


x(n) : Input signal (ampressed voltage or current) 
y(n) : Output signal (electric field , current, etc.) 


We checked that all poles are placed the stability region in z-plane 
by using unknown coefficients a; and bj. 


output y(t) 





output y(t) 


**(1) Z — exp(joT) 
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2.2 Estimation method of Electromagnetic field using a transfer function 





























Input signal ARMA algorithm 
> 
T: 3 
S 9 
S 8 
fam, — 
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Calculated by E es 
FDTD Frequency|Hz| 














Output signal Calculated by time[s] 
ARMA algorithm 
Estimated by 
Transfer function 


E| V/m] 
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3. Simulation result 
3.1 Calculation of the steady state E field 
Analysis Model : Wireless Power Transfer System (WPT System) 


Calculation of the steady state E field when sin wave of 5.5MHz was the input signal 
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F = Resistor =50 Q 
Feed Point 








Calculation point 


WPT System model Skin muscle bone 
skin :¢6=547 [-] o=0.122 [S/m] 
muscle : ¢=282 [-] o=0.594 [S/m] 
bone : ¢ =55.1 [-] o=0.0380 [S/m] |! 
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The calculated data for Transfer function 
Calculation of Transfer function 





H(Z) = YZ _Fly@I] ya): E® [Vim] The E field at each calculation point 
—.X(Z) T[x(t)] xà): Input Voltage [V] oos | | 
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1 3 " 520 v an 
Comparison Result : Transfer function(skin) RS | 


0.020 X Cx "eain. eit bone 
Original FDTD-Real & as skit qusc] 
| ARMA/FDTD-Real Unit : mm 
0.015 Original FDTD-Imag y 
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Calculation point 


-0.015 





frequency [MHz] 


Calculation accuracy : Original FDTD Z ARMA/FDTD 
The number of time step : Original FDTD » ARMA/FDTD 


(240,000steps) — (80,000steps) » 





TAI 


Original FDTD(skin) = 
Original FDTD(muscle) e S —— 
Original FDTD(bone) ~~ iz x id i e [E oua one 
ARMA/FDTD(skin) e. 
ARMA/FDTD(mu Unit: mm 
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Calculation point 
skin muscle bone 





Original FDTD = ARMA/FDTD 
: Original FDTD » ARMA/FDTD 13 
(1,900,000steps)  (80,000steps) 


E TA 


Calculated accuracy : 


The number of time steps 





3.2 Calculation of the Reflection coefficient 
Modell : Loop antenna of WPT system 
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CMS © © Resistor =50 Q 


Tie Feed Point 


Unit : mm 


Model2 : Loop antenna of WPT system near three layered dielectric 


Peed Pori Resistor =50 Q 


Unit : mm 


skin muscle bone 


skin : ¢€=547 [-] o=0.122 [S/m| 
muscle : ¢,=282 [-] o=0.594 [S/m] 
bone  :&,755.1 [-] 070.0380 [S/m] 14 
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"The calculated data for Reflection coefficient 


Calculation of Impedance 


YZ) Fly®] 

H(Z) =— = ——_ 
V) T XG) F[x() 
y(n) : Input Voltage [V] 
x(n) : Input Current [A] 














‘Input Voltage Calculated 
output signal y(t) 
Gaussian pulse 
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The input current at each model 









































Input signal x(t) 
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Comparison Result : Reflection coefficient 


Modell : Loop antenna of WPT system 
Model2 : Loop antenna of WPT system near three layered dielectric 


Modell 


520 


zi = Bu ive 


E * Resistor =50 Q 
~ Feed Point 








Unit : mm 





Reflection coefficient [dB] 
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Calculation accuracy : Original FDTD = ARMA/FDTD 
The number of time steps : Original FDTD >ARMA/FDTD 16 
(1,000,000steps) (100,000steps) 
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4. Conclusion 


Y Our proposed method 
1. ARMA/FDTD method 
2. The estimation method of EM filed using a transfer function 


Y We tested the our proposed methods 
by the calculation of the steady state E field and Reflection coefficient. 
The simulation result 


Calculation accuracy : Original FDTD = Proposed method 
The number of time steps : Original FDTD >> Proposed method 


Our proposed method can reduce calculation time of FDTD method 
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